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Abstract We consider the two dimensional version of a drainage network 
model introduced by Gangopadhyay, Roy and Sarkar, and show that the appro- 
priately rescaled family of its paths converges in distribution to the Brownian 
web. We do so by verifying the convergence criteria proposed by Fontes, Isopi, 
Newman and Ravishankar. 
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1 Introduction and results 

The Two-dimensional Drainage Network Model Let Z be the set of the 

integers and fl = {oj{z), z S Z^) be a family of Bernoulli independent random 
variables with parameter p € (0,1). Let Pp and Ep be the probability and 
expectation induced by those variables in the product space ({0, 1})^^- Consider 
a second family T = ('L'(z), z G Z^) of Bernoulli independent random variables 
with parameter 1/2. Let P1/2 and be the probability and expectation 

induced by these variables in the product space ({0, 1})^^- Let P = Pp x V1/2 
be the product probability induced by the Bernoulli variables in the product 
space {0, 1}^ X {0, 1}^ and let E be the expectation operator with respect to 
this probability. 

For z = (zi, Z2) € 1? 1 say that z' is in the next level if z' — {w, Z2 + 1) for 
some ui G Z. Let h{z) be the closest open vertex to z in the next level with 
respect to the distance induced by the li norm. That it is to say, uj{h{z)) = 1 
^-nd J2i=i \Hz)i- Zi\ = min{^-^^ |a;j-Zi| : uj{x) = 1 and X2 = Z2 + I}. If there 
exist two closest open vertex to z in the next level, the connection will be to its 
left if v{z) = and to its right if v{z) = 1. Now, let h'^{z) = z and iteratively, 
for n > l,/i"(z) = /i(/i"-i(z)). 

Now let G = {V,£) be the random directed graph with vertices V = 7? and 
edges £ = {{u,h{u)) : u This model was proposed in [T] and will be 
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called the Gangopadhyay, Roy and Sarkar (two-dimensional drainage network) 
models or GRS model, for short. 

Main result It follows from the construction that the GRS model may be 
seen as a set of continuous paths, as follows. For any z — (zi, Z2) G we define 
the path X^- — {X^{s), s > Z2} in as the linearly interpolated line composed 
by aU edges {{h^iz),h^+^{z)) : fc e N} of the model, withX^(fc) =h^{z),k€ N, 
where is the fc-th composition of h, meaning the identity. Clearly, X^- is 
a continuous path starting at time Z2- We Let 

X -.^ {X"" : z (1.1) 

which we also call the drainage network , and consider its diffusive rescaling 

Xs {{5xi,5''x2) e : {xx,X2) £ X}, (1.2) 

for 6 G (0, 1]. Our main result below shows that Xs converges in distribution to 
the Brownian web. 

Several authors constructed random processes that formally correspond to 
coalescing one-dimensional Brownian motions starting from every space-time 
point [21 131 m m [6]. In [5], the Brownian web is characterized as a random 
variable taking values in a metric space whose points are compact sets of paths. 
Denote by (T-l,dn) the complete separable metric space where the Brownian 
web is defined. Denote also by T-h the corresponding Borel cj-algebra generated 
by dn- 

The closure in path space of the rescaled drainage network Xs , also denoted 
by Xs, is an (7i, J>^)- valued random variable. 

Theorem 1 The rescaled drainage network Xs converges in distribution to the 
Brownian web as 6 0. 

In order to prove Theorem[l] we will verify the convergence criteria of 5 . To 
describe them, we need the following definition. Given to G R, t > 0, a < b, and 
a (H,J-"-h)- valued random variable V, let r]^{to,t; a,b) be the {0, 1, 2, . . . , cx)}- 
valued random variable giving the number of distinct points in R x {to + t} that 
are touched by paths in V which also touch some point in [a, b] x {to}- 

We can now state the convergence criteria. Let I? be a countable dense set 
of points in R^. 

Theorem 2 (|5j) Suppose that Xi,X2, . . . are {Ti, J--}i)-valued random variables 
with noncrossing paths. If, in addition, the following three conditions are valid, 
the distribution of Xn converges to the distribution of the standard Brownian 
web. 

(Ii) There exist 9f^ G Xn such that for any deterministic yi, . . . ,ym G T>, 
9"^^ , ■ . . , ^fi™ converge in distribution as n ~* 00 to coalescing Brownian 
motions (with unit diffusion constant) starting at yi, . . . , 

(Bi) limsup„^^sup(^ 4^)gR2P(r7^^(to,t;a,a + e) > 2) ^ ase^O+; 
(B2) e"Mimsup„^ocSup(„ t^)gR2 f{-q^^^{to,t;a,a + e) > 3) as e ^ 0+. 
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The present model can be viewed as a discrete space, long range version 
of the two-dimensional Poisson tree introduced in [7] • In [5] , the weak conver- 
gence to the Brownian web of the rescaled collection of paths of that model 
was established. The approach is the same as the one pursued here. There are 
considerable technical differences in the present work, though. The main ones 
are as follows: a) the interaction among the paths in the present model is a 
long range one, differently from the former one; this requires considerably more 
care in the proofs of (Ji) and of an estimate on the tail of the distribution of 
the coalescence time of two given paths, which we use in the proof of {B2) (see 
Theorems [8] and HI , and b) as in ^8], we resort to an FKG type of argument in 
an estimation in the proof of (S2), and for that we need to establish a mono- 
tonicity property of certain conditional probabilities on conditioning paths (see 
Proposition [S]) ; this is also more delicate here than in [S]. 

Another drainage network for which the rescaled paths were shown to con- 
verge in distribution to the Brownian web is the one proposed by Scheidegger 
[9], which is equivalent to coalescing simple random walks. This case was treated 
in [5]. In the general context of drainage networks, our results may be seen as 
part of the scaling theory for those models. See [TU] . 

The GRS model may be seen also as coalescing random walks, in this case 
whose paths are non-crossing and have unbounded dependence among them- 
selves. The other networks mentioned above may also be seen as non-crossing 
coalescing random walks. Coalescing non-simple random walks which are inde- 
pendent before coalescence, and whose paths may thus cross one another, were 
studied in [12] and \TS\ , and shown under suitable conditions to converge to the 
Brownian web when appropriately rescaled. 

2 Proof of Theorem [1] 

This section is entirely devoted to the proof of Thcorcm[T] Let us fix a sequence 
S = Sn = l/\Ai of positive numbers going to zero as rt oo. We want to 
verify conditions Ji , Bi and B2 of Theorem [5] for Xs^ ■ Due to the translation 
invariance of the model, Bi follows from Ii (as will be argued below). So we 
only have to verify B2 and Ii . This will be done in Subsections 12.11 and 12.21 
respectively. As a tool for both this verifications, we start by deriving a bound 
on the time of coalescence of two paths of the drainage model. 

Coalescing random paths Let u,v € be such that u{l) < v{l) and 
m(2) = w(2). Consider AT" and A" and, for t > u(2), define 

Zt^Ztiu,v)^X^-Xl\ (2.1) 

Notice that ^u(2) — v{l) — u{l). 

Remark 3 As argued in ^l] (see 3.1 in the proof of Theorem 2.1), {Zt{u,v) : 
t > u(2)} is a nonnegative martingale in L^. Also, Zt{u,v) —> a.s. as t 00. 



Estimates for the tail of coalescence times We need to control the tail 
of the meeting time of two coalescing random paths starting at the same time 
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at a distance 1 apart. Define 

T := mm{t > : = 0} 



(2.2) 



Theorem 4 There exists a constant C2 > 0, such that P(r > t) < C2/Vt- 

In order to prove Theorem [3] it is enough to consider the case in which 
u — (0,0) and v = (1,0). By Remark|31 Skorohod representation holds: that it 
is to say, there exists a Brownian motion with coefficient of diffusion 1 starting 
from 1 and stopping times = Tq, Ti, ; T2, . . . satifying: 

Zt = B{Tt). (2.3) 

where = Tq, Ti, T2, . . . are such that 

Tt = inf{s>Tt_i :B(s)-B(rt_i) ^ {Ut{B{Tt-i)),Vt{B{Tt-i)))} , (2.4) 

where {{Ut{m), Vt(m)) , t > 1, m G Z} is a famihy of random independent vec- 
tors and for all m e Z, {Ut{ra), Vt(ra)) = {Ut'{m), Vt'{ra)) and {Ut{m), Vt(m)) € 
{{0,0), A} with A = {. .. ,-2,-1} X {1,2,. ..}. Note that since P(Zt > 0) = 1 
we have that 

P {Ui{m) > -m) = 1, V 771. (2.5) 
The following result will be needed later on. 

Lemma 5 For every p < I, there exists a constant ci £ (0,1), such that for 
m > I 

P{iUi{m),Vi{m))^{0,0))<ci. (2.6) 



Proof 

From the Skorohod representation, one readily concludes that the left hand 
side of (|2.6p can be written as P(X^™'°'' — x'f''^'' — m). By conditioning on 
a straightforward computation yields 

as an upper bound for that probability for all 771, where q — 1 — p, and this 
expression is strictly less than 1 for p < 1 . □ 

Proof of Theorem |4] 

Let t' := min{t >0:B{t) = 0}. From equations ([Q]) and (f23|) we get 

P(r > i) = P(r' > Tf)- (2-8) 

Now, for C > 

IP(t' > Tt) < P(t' > Ct) + P{Tt < CO < -7= + P{Tt < Ct), (2.9) 

Vt 
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where c — c(C) S (0,oo), and the second inequality follows by the well known 
behavior of the tail of the distribution of r'. By the Markov inequality, we have 

< C<) = P(e"^^* > e"^'^*) < e^'^* E (e^^^*) . (2.10) 

We now write 

t 

r* = ^5,(z,_i), 

where {Si{k), i > 1, k G are independent random variables. (Note that for 
k G li, {Si{k), i > 1) are not identically distributed.) Then, for A > 



E(e'^^') 



E 



E exp{-A^5,(Z,_i)}exp{-A5t(Zt_i)}|.Ft_i 

\ i=l J 

< E ( exp{-A V 5,(^.-i)} I sup E (e^^^'™)) 



sup . 



where Tt is the cr— algebra generated by {Zq, Zi, . . . , Zt}. 
Hence, 



F(Tt < Ct) < 



supE(e-^5(")) 



(2.11) 



(2.12) 



and by Lemma and Skorohod representation, we get (ommiting m's) 

E(e-^^) < P{{U,V) = {0,0)) 

+ E (e-^^|(C/, V) = (-1, 1)) (1 - P (([/, V) = (0, 0))) 



1 



< 



cosh(V2A) 
1 



P((C/,V) = (0,0)). 



cosh(\/2A) 



cosh(V2A) 



Cl 



cosh(V2A) 



(2.13) 



Indeed, we first notice that 5 is the hitting time of ([/, V) by B, and thus 
the case {U, V) — (—1, 1) is dominated by the cases where {U, V) ^ (0, 0); this 
justifies the first inequality in (j2.13|) . The equality is a well known result (see 
e.g. Theorem 5.7 in llj), and the last inequality follows from Lemma[Sl Since 
Cl < 1 (uniformly in m), we may find Aq > and C^q > such that 



C3 



cosh(v2Ao 



Cl 



cosh(v2Ao) 



< e 



Then, from ([2^2]) 



'(Tt<Cot)<(c3e^''^")* :=c*. 



(2.14) 



(2.15) 



where, from (|2.14p . C4 may be taken < 1. It follows that there exists C5 e (0, 00) 
such that C4 < c^/^/t. Making now C2 = c + C5, the result follows from (|2.9p 
and ^JE\i . □ 
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2.1 Verification of condition B2 

The drainage network model is translation invariant, then we can eliminate 
sup(^_tij>)g22 in {B2) and consider a = tg = 0. 
Then, to prove {B2) it suffices to show that 

e-^limsupP (r]x{0,tN,0,ey/N) > s) ^ 0, if e + . (2.16) 

For j e Zn [0,e//V] let X.j = {X^^'^^{k), keZn [0,tN]} be the [tiVj-step 
trajectory starting in (j, 0) . 

We now introduce the counting variable 77' = | {Xj{tN) 1 < j < n} \ , where 
Xj{tN) ^ X'^^^°'>{ltN\). Then, 

P {r]x{0, tN; 0, eVx) > = P (77' > 3) . (2.17) 

Since 

W > 3} = U]ll{Xj^iitN) < X,{tN) < XnitN)}, 
where n is short for [e\/riJ : have that 

ri-l 

P (??' > 3) < ^ P {Xj^iitN) < Xj{tN) < X„{tN)) 

ri— 1 „ 

= V / P (Xj^iitN) < Xj{tN) < Xn{tN)\Xj = tt) F{Xj = tt) 

n—1 „ 

• P (XjitN) < Xn{tN)\Xj = tt) P{Xj ^ tt), (2.18) 

where lij stands for the state space of Xj . The last equality folows from the fact 
that given X^ = tt, the events {Xj^i{tN) < Xj{tN)} and {Xj{tN) < X„{tN)} 
are independent. 

At this point we want to appeal to the Harris-FKG inequality. So we must 
establish monotonicity properties of the conditional probabilities on the right 

of (ina. 

We begin by introducing a partial order -< on Ilj as follows. Given tti and 
TT2 € Ilj, say that 

TTl ^ TTa ^ TTlil) - 7ri(fc) < TTa (^) - 7r2(fc), (2.19) 

for every £ > k > 0, £, k e 1 D [0,tN]. 

The partial order we have just introduced is an order on the increments of 
the trajectories belonging to events in Ilj. Then, if the increments are increasing, 
the corresponding events will also be increasing. 

Proposition 6 Let 7ri,7r2 G lIj be such that tti -< TT2. Then, for k < j < n 

P [Xk{tN) < Xj{tN)\Xj - TTi] < P [Xk{tN) < Xj{tN)\Xj = ttz] , (2.20) 
P [Xj{tN) < XnitN)\Xj = TTi] > P [Xj{tN) < Xn{tN)\Xj = ttz] . (2.21) 
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Now, using Proposition [SI and since the increments of Xj are independent, we 
may apply Harris-FKG and we find an upper bound for (j2.18p as follows. 

n — l „ 

^ / ¥{X,^,{tN)<X,{tN)\X,^n)¥{X,^7r) 
• / P {Xj(tN) < Xn{tN)\Xj = tt) P{Xj = tt) 

n-l 

= ^ P {X,^i{tN) < Xj{tN)) P {Xj{tN) < X„{tN)) 

ri-l 

< Y.^{X,^i{tN) < X,{tN))F{Xo{tN) < Xn{tN)) 

< nFiXoitN) < XiitN))¥ {Xo{tN) < X^^itN)) 

< eVNF (r > tN) P {t^^n > tN) , 

where t is the coalescing time between the trajectories of X'"'*^) and and 
Te^AT is the coalescing time for X^-^'^^ and X^'^'^K 
From ([rrf)) we get 

e-ip(?7x(0,t7V;0,e%/iV) > 3) < ^/NF {t > tN)P {t^^n > tN) . (2.22) 

From Theorem [51 we have that 

limsupP (t^.at >tN)=V {t,^b > t) , 

where t^^b is the coalescing time between two Brownian motions starting at 
a distance e at time zero. Also, recall that P {t^.b > t) is well known to be 
0(e). On the other hand, by Lemma [H there exists a constant > such that 
P(r > tN) < c'2/VtN. Taking limsup in both sides of (|2?22| we get 

N^oo 

e^MimsupP (r7x(0,iiV;0,e\/iV) > s) = 0(e), 

which shows (|2.16p . 

Proof of Proposition [6] 

It is enough to argue (|2.20p . since by the symmetry of the model we have 

P [Xj{tN) < Xn{tN)\Xj = vr] = P [X2j-n{tN) < Xj{tN)\X, = tt-] , (2.23) 

where 7r~ is the path of Ilj whose increments are the opposite of the respective 
ones of TT. 

The next reduction is that it is enough to consider 7ri,7r2 such that 7ri(s) = 
712(5) for < s < to and 712(5) — tti{s) + 1 for s > to for some to £ [0, tN] D Z. 
From now on, we will be dealing with such pair of paths 7ri and 7r2 . 
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Remark 7 Given that Xj — tt, the distribution of (f2, T) changes as follows. 
Let Iq = (j, 0) and for fc e [1, tN] D Z, let Ik be the interval of integer numbers 
with end points 27r(fc — 1) — 7r(/c) and 7r(fc). (Ik can be a single point.) Then, off 
Ik, the distribution of (fl, T) does not change. In Ik, the variables of 

(f2, T) are independent of the ones off U^^q Ik, and they are also independent 
among themselves, except the pairs (ct;(27r(fc — 1) — 7r(fc), k), v{'K{k — 1), fc — 1)) 
for which \Ik\ > 1, whose distribution is described below. First, for each k G 
[IjtA^] n Z, u}{-K{k),k) — 1. // \Ik\ > 1, then w(zi,fc) = for z\ in Ik minus 
its endpoints, and u;(27r(fc — 1) — 7r(fc),fc) is a Bernoulli random variable with 
parameter p' — p/ {2 — p); given that uj{2n{k — 1) — 7r(fc), k) = 0, then v{'K{k — 
1), fc — 1) is Bernoulli with parameter 1/2; otherwise it is Q or 1 depending on 
whether ■K{k) — n{k — 1) > or < 0. 

A further point is that given that Xj = n, the event {Xk(tN) < Xj{tN)} 
depends on i7 only to the left of and including tt, and on T only strictly to the 
left of TT . 

In order to prove ()2.20p . we couple the distributions of (11, T) given vri and 
TT2 conveniently (in the relevant region of I?). Let us denote by {uJi{z),Vi{z)), 
z S Z^, the general term of (fi, T) given tt^, i — 1,2, respectively. Then 

{uj2{z),V2{z)) = {uji{z),vi{z)) for < Z2 < to, and {uj2{z),V2{z)) = (wi((zi - 
1, Z2), vi{zi — 1, 22)) for < Z2 > to- The coupling on Z2 = to will be different 
in different cases. 

Case 1: 7ri(to) < 7ri(to — 1). 

In this case, we make (a;2(-zi, to), ^2(^1, to)) = {u}i{zi — l,to),ui(zi — l,to)) 
for zi < 7r2(to), and (^2(2^1, to), ^2(21, to)) = {oJi{zi + 1, to), 171(21 + l,to)) for 
zi > 7r2(to). 

One readily checks in this case that for every realization (wi, ui)(-) of (51, T) 
given Xj = tti, {lo2, W2)(0 given by the above coupling is a realization of (il, T) 
given Xj — 112, and that 

{{uJi,vi){-) : Xk{tN) < X,{tN)} C {{co2,V2){-) : Xk{tN) < Xj{tN)}. (2.24) 

This then establishes (|2.20p in this case. 

Case 2: 7ri(to) > 7ri(to - 1). 

We start with a coupling which is the symmetric of the one above. It does 
not quite give us (I2.24|) . but something nevertheless useful. 

Auxiliary coupling: 
Given (tj2,U2)(-), we make (wi(zi, to), t;i(2i, to)) = (0)2(21 + 1, to), ^2(2:1 + 1, to)) 
for zi > 7ri(to), and (wi(2i, to), ui(zi, to)) = (^2(21 - 1, to), t;2(^;i - l,to)) for 
zi < 7ri(to). One can readily check that this provides a coupling of (il, T) 
given Xj — tti and (11, T) given Xj — 7:2. It is possible to find realizations of 
(w2,W2)(') for which with this coupling (|2.20p does not hold. But the following 
can be readily checked. 

{(wi,wi)(-) : Xkito) < Xjito)} C {iL02,V2)i-) : Xk{to) < Xj(to)}. (2.25) 

This will be used in the coupling discussed next. 

We begin by pointing out that we could not indeed find a direct coupling 
for which (|2.24p holds, so we take an indirect route. We start by considering 
Dq = Xj{tQ) — Xk{to). In order to insure the existence of a coupling of (wi, ui)(-) 
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and {lu2, W2)(-) satisfying (|2.24p . it is enough to show that the distribution of Do 
given Xj = tt2 dominates that of Dq given Xj — tti. Indeed, if this holds, then 
we can find a couphng of the two distributions such that the domination takes 
place almost surely. The respective random variables depend of course only on 
the history up to to, and thus we are free to choose any coupling of {uji,vi){-) 
and {ui2,V2){-) above to, and choosing as in Case 1, we readily get (|2.24p . 

In order to establish the above mentioned domination, we first define A = 
2(71-1 (to) -TTi (to -1)) and Dq = 7ri(to) - A - Xfc(to). Denote by A, « = 1,2, the 
random variable whose distribution equals that of the conditional distribution 
of Dq given that Xj — tt^, i — 1,2, respectively. Then the above mentioned 
domination is equivalent to 

F{D2>k-l)>F{Di>k),k>0. (2.26) 

We will argue ()2.26p using different couplings for different cases. The case 
A; = of (|2.26p is established as follows. If 7ri(to) = 7ri(to — 1), there is nothing 
to prove. If 7ri(to) > 7ri(to — 1), then _Di > is equivalent to the left hand side 
of (|2.25p . and > ^ 1 is equivalent to D2 > 0, which is in turn equivalent to 
the right hand side of ()2.25p . So the auxiliary coupling can be used in this case. 
It can also be used if fc = 1. Indeed, in this case, D2 > is equivalent to the 
right hand side of (|2.25p and Z?i > 1 is either equivalent to the left hand side 
of (|2.25p (if 7ri(to) = ■7ri(to — 1)), or is in any case contained in it. 

For k > 2, we will resort to another coupling, described as follows. 

1. For zi > 7ri(to), {uj2{zi,ta),V2{zi,to)) = {uji{zi - l,to),vi{zi - l,to)); 

2. for 7ri(to) - 2(7ri(to) - 7ri(to - 1)) < 21 < 7ri(to), ^2(^1, to) = and 
^2(21, to) are i.i.d. BernouUis with parameter 1/2 independent of all else; 

3. for zi = 7ri(to) - 2(7ri(to) - 7ri(to - 1)) - 1, (^2(^1, to), ^2(21, to)) = 
{^LUi{zi,to),v), where ^ is a Bernoulli with parameter 1/(2 ~ p) inde- 
pendent of all else, and v depends only on ^ll!i{zi, to); if the latter random 
variable vanishes, then v is Bernoulli with parameter 1/2 independent of 
all else; otherwise, v = 1; 

4. and for zi < 7ri(to) - 2(7ri(to) - 7ri(to - 1)) - 1, (1^2(21, to), ^2(21, to)) = 
{uJi{zi,to),vi{zi,tQ)). 

One readily checks that for every realization {uji,vi){-) of (il, T) given Xj = 
TTi, {lu2, V2){-) given by the above coupling is a realization of {Q, T) given Xj — 
7r2, and that 

{{L^i,vi)i-) : Di>k}ci {iL02,V2)i-) : D2>k} (2.27) 
in the remainder cases. □ 

2.2 Weak convergence to coalescing Brownian motions 

The main purpose of this part of the paper is to prove condition II of Theorem 
m Also, we prove condition Bl, which is in fact a consequence of condition II 
for this model. 

Denote a single trajectory belonging to the drainage network model and starting 
from a e Z at time j £ Z hy X"^^ = {X°'^\l),l > j}. Then, 
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Theorem 8 Let {yo, sq), (yi, si), . . . , (yk, Sk) be fc + 1 different points in such 
that So < si < ■ ■ ■ < Sk and if Si_i ~ Si for some i,i — 1, . . . ,k, then yi-i < yi- 
IfZ^n^ = {Zli\t) := n-i/2x«([nij) := n-i/^XL^^-v^J^L-'^J ([ntj), i > sj; then 

{Z^\i = 0, . . . , fc} ^ {W^W, z = 0, . . . , fc}, (2.28) 

as n ^ oo, where the W^^"^ are k + 1 coalescing Brownian motions with positive 
diffusion coefficient a , starting at {{yo, So), • . ■ , {yk, **;)}• 

Remark 9 Here, =^ stands for weak convergence in 11*''+^ (the product of k + 1 
copies o/n, the path space (see JS^)- 

Remark 10 is the variance of X'^'^{1), which one readily computes from its 
distribution, which is in turn straightforward to obtain, as ^-i^^^^yi ■ 

Proof of Theorem M 

The proof is divided in three parts. The arguments are somewhat standard, 
so we will be sketchy at a few points. 

Part I 

When fc = 0, the result follows from Donsker' invariance principle for single 
paths since X"-'^ is a random walk with variance . 

Part II 

In this part, we consider the case where fc > 1 and so — ■ ■ ■ ~ Sk (and 
yi-i < yi for every i,i — 1, . . . , fc). The general case will be treated on Part 
III. Since the increments are stationary and spatially homogeneous we will take 

so^yo^ 0. 

a) We first treat the case fc = 1 . Without loss of generality, we suppose that 
we are working in a probability space in which 

ZW^J!^W^^\ (2.29) 

n — >co 

The following proposition, besides implying the result in the case fc = 1, is also 
the building block of the argument for the other cases of this part. 

Proposition 11 The conditional distribution of ^ given converges al- 
most surely to the distribution of W'^'^^ given 

Proof 

The strategy is to approximate the paths Zn'^ and zl^^ before the time when 
they first come close (at a distance of order n", a < 1/2) by independent paths; 
then to show that after that time they meet quickly. Let us write 

Xi^\l)^[y,V^\+Y.Sh"\ (2.30) 

where for each i — 0,1 and n > 1, sl^'^\h > 1, are i.i.d. random variables 
distributed as X"'°(l). Next, we introduce independent copies of sl^'^\i — 
0,1, n > 1. For each n > 1, let {sl^'"'\h > l,i — 0,1} be i.i.d. random 
variables such that s[^'"^ ^ X°'"(l) and define random variables sj^'"^ by 

I 5^ , otherwise. 
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where a is a positive number to be specified later. Moreover, for Z > 1, let 
X«(0 = Ly.x/^J +EL=i^I''"^- Define 



Tn = 

and 

Tn = 

and let 



niin{Z > 1 : Xj^^\l) - X^°\l) < 3n"} (2.32) 
niin{Z > 1 : X^^\l) - X^^^l) < 3n"}, (2.33) 



X^^ {Tn) + J2h> fr, ^h'"' ' Otherwise, 



where 



^ 1 Sj^'"-\ ifi = l. ^ ^ 



ifi = l 

-,(i,n) 



Prom the exponential tail of the distribution of 5^ ' ' (which is independent 
of i) we readily get the following result. 

Lemma 12 Given (3 <oo, the complement of the event 

An :={4*'"^ = 5(*'"\i = l,...,n^i = 0,l}, 

has probability going to zero superpolinomially fast as n ^ oo. 

Corollary 13 The event := {xi,'\l) = X^\l);l = 1, . . . , n^; i = 0,1}, is 
such that P[i?fJ converges to superpolinomially fast as n ^ oo. 

From the construction, the properties stated in the following result are read- 
ily verified. 

Lemma 14 X^n ^ {xi°\l), I > 1} and xi^^ := {^^^'(0, I > 1} are indepen- 
dent processes and Xn^ ^ Xn^ for each fixed i = 0,1. 

Corollary 15 Let Z^\t) := n'^/"^ X^\lnt\), t > 0. Then, 

1. Zn^ W^^^ almost surely as oo; 

2. {zi^\ziP) ^ (W(o), M^(i)), where ~ W^, and and W'^^^ are 
independent; 

3. the conditional distribution of {zi^\t), t < ^} given ZrP converges al- 
most surely to that o/{H^^°^(t), t < t} given W^-'-\ where r — inf{t > : 

w(-°\t) = w(-^\t)}. 

Corollary 16 The conditional distribution of {Zn\l),l < Tn} given Zn ^ con- 
verges almost surely to that of {W^^\t), t <t} given W^^h 

Lemma 17 For a < 1/2, there exists 7 = 7(a) < 1 such that, as n ^ 00, 

P[X(0)(0 ^ {l),l = Tn,...,Tn+ n-^] 0. (2.36) 

Corollary 18 {zL°\t),t>^} converges in probability to {W''^\t), t> t}. 
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Then, PropositionfTDfoUows directly from CoroUariesfTBlandfTSl and the Markov 
property. □ 

Proof of Corollary 1151 

1. Let /3 > 7/2 be fixed. We claim that there exists a finite constant C and 
9 > \ such that 

P[t„ > nf^] < Cjn^. (2.37) 

To argue this we first notice that := X^^^ - xf^ is a random walk on 

the integers starting from [yiv^J ^'^^ with increments distributed as F("^(1) = 
^(i,n)_^(o,„)^ One may readily notice then that f„ is dominated by the sum of 
Lz/iv^J copies of f„, the hitting time of the negative integers by a random walk 
with increments distributed as y'^"^(l), started at the origin. We then have that 

P(f„ > n^) < ci V^P(f„ > C2 r/"i/2) < C7/n(''-3/2)/2, (2.38) 

where the first bound is straightforward from the above domination, and the 
last one follows from a standard fact on the tail of the distribution of hitting 
times for one dimensional random walks with finite second moment, like Y^"'\ 
The claim is thus justified. 

We can then use ([T^ . CoroUary [HI and (fnTj) . together with Borel- 
Cantelli, to conclude. 

2. Immediate from Lemma [17] and Donsker's theorem. 

3. We take a version of Zn\ say Zn\ such that 

(Z(°),ZW)^(1^(°),1F«), (2.39) 

with {W^°\W^'^'>) ^ (W^(°),M^(i)), and claim that f„ ^ f almost surely as 
n ^ oo, where f„ = min{i > : Z^^\t) - zi^\t) < 3n"-i/2} and f = inf{t > 
: VF(o)(i) = 1^(1) (i)}. Indeed, the almost sure convergence of r„ follows from 
the property, that with probability one, two independent Brownian trajectories 
will instantaneously cross after touching, i.e. 

P[Ve > 0,W^°\t) > W'-^^t) for some t<f + e] = l, 

and the following deterministic result (whose proof is an exercise). The conver- 
gence in distribution in the original space follows. □ 

Lemma 19 Let fn, gm f, 9 ■ R be continuous functions such that /(O) < 

g{0), T = inf{t > : f(t) — g{t)} is finite and has the property that for every 
S > there exists t G [T,T + S] with f{t) > g{t). Suppose in addition that 
lim„^oosupo<i<j.+i |/„(t)-/(t)| =0, anrf lim„^oo supo<t<-r+i \gn{t) - g{t)\ = 
0, and let T„ = inf{i > : gn{t) — fn(t) < e„}, where (e„) is a given sequence 
of numbers vanishing as n —> oo. Then, lim„^oo Tn = T . 

Proof of Corollary 1161 

Immediate from CoroUarv fT51 and 3. of Corollary [131 □ 

Proof of Lemma 1171 

Let cr„ — ini{l > Tn '■ X^\l) — Xn \l)} ~Tn- It can be readily checked that 
(Tn is dominated by an, the hitting time of by Zt{{0, 0), ([3ri"J , 0)) (see (|2.ip ). 
Proceeding as in the proof of Theorem [H we find 

P(a„ > n^) = P(a:, > Tn-,) < P« > Cn^) + P(r„. < Cn''), (2.40) 
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where ct^ is the hitting time of by a Brownian motion started at [3n"J 
(see (|2.8p and (|2.9p ). Choosing now 7 in (2a, 1), one readily checks that the 
first term on the right of (|2.40p vanishes as n ^ cxd. That so does the second 
one follows by ()2.15p . and the proof is complete. □ 

Proof of Corollary [18] 

By (I2.29P and Lemma [T71 it is enough to show that for every e > 0, as n ^ 

P( max {Xi^^l + r„) - + r„)) > eV^) ^ 0. (2.41) 

0</<nT 

By the Markov property and an elementary domination argument, the left 
hand side of (|2.4ip is bounded above by 

P( max Z,'"^ > e^/H,), (2.42) 

0<i<n-' ' ' 

where Z^"'' = Zi{(Q,0), ([3n"J,0)) (see ^^). Now by Doob's inequality, (|2^ 
is bounded above by constant times E(Z,^"'')/\/n < 3n"^^/^ ^ as n ^ 00. □ 

b) We now consider the case A; > 1. It suffices to prove that 

E[/o(z(o)) . . . A(ZW)] > E[/o(W-(°)) . . . /fe(W-«)] (2.43) 

n — >oo 

for any fo, ■ ■ ■ ,fk G Cb(n, Hi), the space of all real- valued and bounded functions 
defined on 11 . 

The left hand side of l|2.43p is equivalent to 

E{A(Z,W)E[/o(Z("))/2(Z(2)) . . . /,(ZW)|Z,W]}. (2.44) 

It follows from Lemma [50] below that the conditional expectation equals 

E[/„(Z("))|ZW] E[/2(Z(2)) . . . /,(ZW)|ZW]. 

Since Zn'' -*■ M^^^-* by the inductive hypothesis the conditional expecta- 
tions above converge, almost surely, to E[fo{W'^°^)\W'^^^ and E[/2(PF(2)) . . . 
/j.(VF('=))|M^(i)], respectively. On the other hand, by the dominated conver- 
gence theorem we get the convergence of (I2.44p to 

E{/i(W^(i))E[/o(W^("))|M^(i)] e[/2(w^(^)) . . . fkiW^'^W^'^]}, 

which equals 

E{/i(W^(i))E[/o(W^(''))/2(W^(')) . . . fkiW^'^^W'-^ (2.45) 
Since (|2.45p equals the right hand side of (|2.43p . we have the result. □ 

Lemma 20 For each n > 1, given Zn \ Zn^ is independent of {Zn \ ■ ■ ■ , zl^^). 
Proof 

Zn^ depends only on the variables of (il, T) to the left of and including 
z'iP , and {Z^\ . . . , zH^'') depends only on the variables of (fi, T) to the right 
of and including Zn'' ■ The result follows from the product structure of the 
conditional distribution of (fi, T) given Zn'' (see Remark[7]), and the facts that, 
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given Zn \ the variables of along Zn'' are fixed, and when Z^'' depends on 
the variable of T in a particular position of ZrP (that happens when, and only 
when, the increment of Zn"^ from the given position is strictly positive), then 
{Z^\ . . . , Zri'^'') does not depend on that particular variable, and vice- versa. □ 

Part III 

We now take on the general case and proceed by induction on k > 0. The 
case k = was argued in Part I above. Suppose now that for some fc > 1, the 
result holds up to fc — 1, and let (yo, sq), . . . , (yk, Sk) be as in the statement of 
Theorem [H If the condition on (s^) of Part II holds, then there is nothing to 
prove. Otherwise, there exists i € {I, . . . , k} such that Si > Si_i. Let iq be the 
maximum between such i's. By the induction hypothesis, 

s,<t< s,„}, 0<i<io)^ {{W^"\t), s,<t< s,„}, 0<i< zo), 

(2.46) 

and we can suppose that such convergence occurs almost surely. By the Markov 
property and Part II (strengthened to accomodate the case where the y^'s may 
depend on n and converge to, say, jji as n —> oo; this requires minor changes in 
the argument above), we have that, given ()2.46p . 

({Z«(i), t>s,,},0<i<k + l)^ ({VK« it),t>s,,},0<i<k + 1). (2.47) 

Then, Theorem [S] follows from (jTIS)) . ([TiT)) and the Markov property. □ 

Verification of Bi 

We note that rix„ (0, t; 0, e) > 2 if and only if the two trajectories starting 
at time from the border of the interval [0, e] have not met up to time t. By 
Theorem |S1 the rescaled trajectories of a finite collection of paths belonging to 
the drainage network model converge to those of coalescing Brownian motion. 
Then, it is straightforward to get, for the drainage network model, that 

limsupP(r/A'„(0,t;0,e) > 2) = 2$(e/\/2l) - 1, 

n — *C30 

where $(•) is the standard normal distribution function. This together with the 
translation invariance of the model concludes the proof. 
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